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We consider Hamiltonian PDEs that can be split into a linear unbounded 
operator and a regular non linear part. We consider abstract splitting methods 
associated with this decomposition where no discretization in space is made. We 
prove a normal form result for the corresponding discrete flow under generic non 
resonance conditions on the frequencies of the linear operator and on the step size. 
This result implies the conservation of the regularity of the numerical solution 
associated with the splitting method over arbitrary long time, provided the initial 
data is small enough. This result holds for numerical schemes controlling the 
round-off error at each step to avoid possible high frequency energy drift. We 
apply this results to nonlinear Schrodinger equations as well as the nonlinear 
wave equation. 
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1 Introduction 

In this work, we consider a class of Hamiltonian partial differential equations 
whose Hamiltonian functions H = Hq+P can be divided into a linear unbounded 
operator Hq with discrete spectrum and a non linear function P having a zero 
of order at least 3 at the origin of the phase space. Typical examples are given 
by the non linear wave equation or the non linear Schrodinger equation on the 
torus. 

Amongst all the numerical schemes that can be applied to these Hamiltonian 
PDEs, splitting methods entail many advantages, as they provide symplectic and 
explicit schemes, and can be easily implemented using fast Fourier transform if 
the spectrum of Hq expresses easily in Fourier basis. Generally speaking, a 
splitting schemes is based on the approximation 

^'h^VHo°^P (1-1) 

for small time h, and where (/?^ denotes the exact flow of the Hamiltonian PDE 
associated with the Hamiltonian K. 

The understanding of the long-time behavior of splitting methods for Hamil- 
tonian PDEs is a fundamental ongoing challenge in the field of geometric inte- 
gration, as the classical arguments of backward error analysis (see for instance 
[16]) do not applied in this situation where the frequencies of the system are 
arbitrary large, and where resonances phenomenon are known to occur for some 
values of the step size. Recently, many progresses have been made in this di- 
rection. A first result using normal form techniques was given by Dujardin & 
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Faou in [6] for the case of the hnear Schrodinger equation with small potential. 
Concerning the non linear case, results exists by COHEN, Hairer & LuBiCH, 
see [8, 9], for the wave equation and Gauckler & Lubich, see [10, 11], for 
the nonlinear Schrodinger equation using the technique of modulated Fourier 
expansion. However to be valid these results use non-resonance conditions that 
are generically satisfied only under CFL conditions linking the step-size h and 
highest frequencies of the space discretization of the Hamiltonian PDF. 

In this paper, we use a normal form techniques to prove the long-time preser- 
vation of the regularity of the initial solution under generic non resonances condi- 
tions valid for a large set of equations and a large set of time step, independently 
of the space discretization parameter. 

Normal form techniques have proven to be one of the most important tool 
for the understanding of the long time behaviour of Hamiltonian PDF (see [1, 
4, 12, 2, 3, 13]). Roughly speaking, the dynamical consequences of such results 
are the following: starting with a small initial value of size e in a Sobolev space 
i?^, then the solution remains small in the same norm over long time, namely 
for time t < Cr£~^ for arbitrary r (with a constant Cr depending on r). Such 
results hold under generic non resonance conditions on the frequencies of the 
underlying linear operator Hq associated with the Hamiltonian PDF, that are 
valid in a wide number of situations (nonlinear Schrodinger equation on a torus 
of dimension d or with Dirichlet boundary conditions, nonlinear wave equation 
with periodic or Dirichlet conditions in dimension 1, Klein Gordon equation on 
spheres or Zoll manifolds.). 

In this paper, we mainly show that the same kind of results hold true for 
numerical solutions associated with the abstract splitting method (1.1) under 
some further restrictions specifically induced by the time discretization. 

This work is the second of a series of two. 

In the first part [7], we consider full discretizations of the Hamiltonian PDF, 
with a spectral discretization parameter K. We show that under the hypothesis 
K < £~°' for some constant a depending on the precision degree r, the same 
conclusion as in the continuous case can be drawn. Though concerning the 
discretization of a Hamiltonian PDF, the method used in [7] is by essence a 
finite dimensional Birkhoff normal form result using techniques that are rather 
classic in the dynamical system world. 

In some sense, the present paper studies the case where K > e~'^ by consider- 
ing the splitting method where no discretization in space is made (i.e. K = -|-oo). 
The techniques used involve the abstract framework developed in [4, 12, 2]. How- 
ever, instead of being valid for the (exact) abstract splitting (1.1), we have to 
consider rounded splitting methods of the form 

H^,, o o (1.2) 

where H^^^ puts to zero all the frequencies whose weighted energy 
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in the Sobolev space is smaller than a given threshold ry^. Hence, for small rj, 
(1.2) is very close to the exact splitting method (1.1). The good news is that this 
threshold can be taken of the order e^, making this projection 11^ ^ very close to 
the identity, and in any case producing an error that is far beyond the round-off 
error in a computer simulation (particularly for large s). 
Our main result is given by Theorem 4.2. 

2 Description of the method 

Before going on into the precise statements and proofs, we would like to give 
tentative explanations of the modifications observed in comparison with the con- 
tinuous case. 

The method used in [4] to prove the long-time conservation of Sobolev norms 
for small data is to start from a Hamiltonian H = Hq+P depending on an infinite 
number of variable (Cj,%)) j,k £ N, and for a fixed number r, to construct a 
hamiltonian transformation r close to the identity, and such that in the new 
variable, the Hamiltonian can be written 

Ho + Z + R (2.1) 

where Z is a real Hamiltonian depending only on the action Ij = ^jr]j and R a 
real Hamiltonian having a zero of order r. 

The key for this construction is an induction process where, at each step, the 
solution of an homological equation of the form 

{H^,X] + Z = G (2.2) 

where G is a given homogeneous polynomial of order n, and where Z depending 
only on the action and x ^'^^ unknown. Assume that G is of the form 

G = Gjk ^ji • • • CjpVki ■■■Vkg 

where Gju is a coefficients, j = (ji,...,jp) G and k = {ki, . . . , kg) G 
Then it is easy to see that the equation (2.2) can be written 

^{j, k)xjk + Zjk = Gjk (2.3) 

where 

k) =ujj^-\ h ujjp - Wfci ujj^ 

and where Zjk and Gjk are unknown coefficients. 

It is clear that for j = k (up to a permutation), we have ^{j, k) = which im- 
poses Zjk = Gjk- When j ^ k (taking into account the permutation), the solu- 
tion of (2.3) relies on a non resonance conditions on the small divisors ^{j, k)~^. 
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In [4], Bambusi & Grebert use a non resonance condition of the form 

Vj/fc, \n{j,k)\>jfi{j,kr^ (2.4) 

where k) denotes the third largest integer amongst . . . , |A;g|. They more- 
over show that such a condition is guaranteed in a large number of situations 
(see [4], [12] or [2] for precise results). 

Considering now the splitting method (/j^^ o ipp, we see that we cannot work 
directly at the level of the Hamiltonian. To avoid this difficulty, we embed the 
splitting into the family of applications 

[0, 1] 9 A ^ ip'h, o V9^p 

and we derive this expression with respect to A to work in the tangent space, 
where it is much more easy to identify real Hamiltonian than unitary flows. 

This explains why we deal here with time-dependent Hamiltonian. Note that 
we do not expand the operator 99^ in powers of h, as this would yields positive 
powers of the unbounded operator Hq appearing in the series. Unless a CFL 
condition is employed, this methods do not give the desired results (and do not 
explain the resonance effects observed for some specific values of h, see [7]). 

Now, instead of (2.2), the Homological equation appearing for the splitting 
methods is given in a discrete form 

Xoy^Ho-X + Z = G. (2.5) 

In terms of coefficients, this equations yields 

The main difference with (2.3) is that we have to avoid not only the indices 
{j,k) so that i}{j,k) = 0, but all of those for which hil.{j,k) = 2rmr for some 
(unbounded) integer m. 

In the case of a fully discretized system for which Vz-P = for \j\ > K, 
then under the CFL-like condition of the form hK^ < C for some constant C 
depending on r, then we have \hQ.{j, k)\ < vr, and hence 

|gih,t7y,fc) _ ^1 > h-^^^{^, fc)-" (2.6) 

(2.6) is then a consequence of (2.4). Under this assumption, we can apply the 
same techniques used in [4] and draw the same conclusions. This is the kind of 
assumption made in [8] and [11]. 

The problem with (2.6) is that it is non generic in h outside the CFL regime. 
For example, in the case of the Schrodinger equation, the frequencies of the 
operator Hq are such that ujj ~ j^. Hence, for large iV, if (ji, . . . , jp, fci, . . . , kq) 
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is such that ji = N + 1, ki = N and ah the other are of order 1 (N is large here), 
we have fc) ~ (A^ + 1)^ - A^^ ~ 2N. Hence, 

|gih.n(j,fc) _ ]_| ^ |g2ihJV _ -,^1 

cannot be assumed to be greater than hjfi^j ,k)~"' ~ h. Note that a generic 
hypothesis on h would be here that this small divisor is greater than h'yN~^ for 
some constants 7 and a. This means that we cannot control the small divisors 
|gj/if2y,fc) _ -|^| associated with the splitting scheme by the third largest integer in 
the multi index, but by the largest. 

Using a generic condition on /i < /iq, we can prove a normal form result and 
show that the flow is conjugated to the flow of a hamiltonian vector field of the 
form (2.1), but where Z now contains terms depending only on the actions, and 
supplementary terms containing at least two large indices. Here, large means 
greater than e"*^ where a depends on r. 

In the case of a full discretization of the Hamiltonian PDE with a spectral 
discretization parameter K, we thus see that if X < e~'' then the normal form 
term Z actually depends only on the actions, as the high frequencies greater that 
e~'^ are not present. This is essentially the result of [7]. 

In the case where K > e"'^ , the normal form result that we obtain can be 
interpreted as follows: the non conservation of the actions can only come from two 
high modes (of order greater than e~^) interacting together and contaminating 
the whole spectrum. The role of the projection operator H^^s is to destroy these 
high modes at each step but only when these high modes have an energy greater 
than Tj (cf. (3.13)). As we can take rj = e"^ , the error induced is very small, 
and in particular, far beyond the round-off error in the numerical simulation. 
Note that the complete understanding of the numerical phenomenon, and in 
particular the possible interaction (or not) between two high modes would in 
principle require the introduction of round-off effects by adding stochastic terms 
in high frequencies. This is clearly out of the scope of this paper. 

3 Setting of the problem 

3.1 Abstract Hamiltonian formalism 

We denote = Z"' or N'^ (depending on the concrete application) for some d > 1. 
For a = (ai, . . . , ad) G M, we set 



1 9 / 9 9 \ 

a\ = max (1, a^^ + • • • + a^) . 



We consider the set of variables {(,a, Vb) £ C-'^ x C-'^ equipped with the symplectic 
structure 
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We define the set Z = N ^ For j = (a, 6) G we define \ j\ = \a\ and we 

denote by j the index (a, —5). 

We will identify a couple G x with {zj)j^z G via the formula 



{a, 5) G Z 



if 5=1, 
r/a if 5 = -1, 



By a slight abuse of notation, we often write z = rf) to denote such an element. 

For a given real number s > 0, we consider the Hilbert space Vg = C) 
made of elements z G such that 



and equipped with the symplectic form (3.1). 

Let Z// be a an open set of Vg- For a function F of C^{U,C), we define its 
gradient by 

'dF\ 



VF z) 



where by definition, we set for j = {a, 6) £ M x {±1}, 

0-, ^ if i = 

Let -^(2;) be a function defined on U. If is smooth enough, we can associate 
with this function the Hamiltonian vector field Xf{{z) defined by 

Xh{z) = JVH{z) 

where J is the symplectic operator on Vg induced by the symplectic form (3.1). 
For two functions F and G, the Poisson Bracket is defined as 

We say that 2; G "Ps is real when z-j = 'Zj for any j G Z. In this case, 2; = (^, ^) 
for some i G C^. Further we say that a Hamiltonian function H is real if H{z) 
is real for all real z. 

Definition 3.1 Let s > 0, and letU be a neighborhood of the origin in Vg- We 
denote by 7i^{V() the space of real Hamiltonian H satisfying 

HeC°°iU,C), and Xh & C°^{U,Vg). 
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With a given function H G 7i^{V(), we associate the Hamihonian system 

i = JVH{z) 



which can be written 



OH 
orja 



(3.2) 



In this situation, we define the flow ip\j{z) associated with the previous system 
(for times t>0 depending on z &U). Note that if z = (^,^) and using the fact 
that H is real valued, the flow (^^^r/*) = ^p^^{z) satisfies for all time where it is 
defined the relation ^* = ff , where ^* is solution of the equation 



dH — 

OVa 

In this situation, introducing the real variables Pa and qa such that 



(3.3) 



1 



1 



= —7={Pa + iqa) and = —7={Pa - iqa 



V2 

the system (3.3) is equivalent to the system 



Pa 
Qa 



-K—{q,P) aeM 
oqa 



dH 

dpa 



{q,p), aeJ\f. 



Note that the flow r* = y?* of a real hamiltonian x defines a symplectic map, 



where H{q,p) = H{(,,(^) 
Note that the flow t' 
i.e. satisfies for all time t and all point z where it is defined 

(L',r*)f J(L),r*), = J 



(3.4) 



where denotes the derivative with respect to the initial conditions. 
The following result is classic: 



Lemma 3.2 Let U and W he two domains ofVg, and let t = tp^ e C°°{U,W) 
be the flow of the real hamiltonian x- Then for K G 'H'^(yV), we have 

yz^U XkoAz) = {D,t{z))-^Xk{t{z)). 

Moreover, if K is a real hamiltonian, K or is a real hamiltonian. 
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3.2 Function spaces 

We describe now the hypothesis needed on the Hamiltonian H. 

Let £ > 3 be a given integer. For j = (ji, . . . ,jr) G Z^, we define fi{j) as the 
third largest integer between | ji | , . . . , | I • Then we set S{j) = \ — \ 
where and \ denote the largest and the second largest integer between 

Let £ > 3. We consider j = {ji, ■ ■ ■ , j^) S Z^, and we set for all i = 1, . . .p 
ji = {ai,6i) where Ui £ Af and 6i G {±1}- We define the moment of the 

multi-index j by 

M{j) = ai5i + ■ ■ ■ + aeSe. (3.5) 
We then define the set of indices with zero moment 

Ii = {j = {ji,...,ji)£Z', with Mij) = 0}. (3.6) 

In the following, for j = (ji, . . . , j^) G 2^, we use the notation 



Moreover, for j = (ji, . . . ,je) £ le with ji = (aj, 6i) £ M x {±1} for i = 1 
we set 

3 = Gi,---,Je) with J. = (oi, -(5i), i = !,...,£. 
We recall the following definition from [12]. 

Definition 3.3 Let k > 3, M > and v G [0,+oo), and let 

k 

We say that Q G 7'*^''^ {f there exist a constant C depending on M such that 



y£ = 3,...,k, yjele, \Qj\ < C^]:/ . (3.7) 



SUV 



Note that Q is a real hamiltonian if and only if 

yi = 3,...,k, yjeli, Qj = Qj. (3.8) 

We have that T,^'^'" £ Ti" for s > u + 1/2 (see [12]). The best constant in the 
inequality (3.7) defines a norm [QI^a/,^ for which '71*^''^ is a Banach space. We 
set 

Men 
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Definition 3.4 A function P is in the class T if 

• P is a real hamiltonian and exhibits a zero of order at least 3 at the origin. 

• There exists sq > such that for any s > sq, P G 7i'^{U) for some neigh- 
borhood U of the origin in Vg ■ 

• For all k > 1, there exists > such that the Taylor expansion of degree k 
of P around the origin belongs to T^''^ ■ 

With previous notations, we consider in the following Hamiltonian functions of 
the form 

H{z) = Ho{z) + P{z) = uJaUz) + P{z), (3.9) 

aeA/" 

where for all a G A^, 

Ia{z) = iaVa 

are the actions associated with a £ M and where G M are frequencies satisfying 

VaGAA, \uJa\<C\ar (3.10) 

for some constants C > and m > 0. The Hamiltonian system (3.2) can hence 
be written 

dP 

= -iiOa^a-i-K {S,,V) CL^M 

Qp'' (3.11) 
r]a = ii0ar]a + iT7^{(,,v) a£M. 

3.3 Rounded splitting methods 

When considering the numerical simulation of such hamiltonian system, many 
methods can be interpreted as splitting methods associated with the decompo- 
sition (3.9). This means that for small step size h, we approximate the flow (p^ 
by the composed flow 

For a given time t, and a small step size h with t = nh, the approximation of 99^ 
is then written 

^^HC^I^^'h.o^'^y . (3.12) 

We give examples of such schemes in the next section. 

In order to control the possible numerical instabilities due to the interaction 
of high frequencies, we introduce the following projection operator: Let rj > 
and s be given, we define 
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by the formula 

/ N { Zj if \3\"\zj\ <V ^ ^ 

Vj G Z, (Ur,,sz) = { (3.13) 

[ if \j\'\zj\ > rj. 

The goal of this paper is the studying of the long-time behavior of rounded 
splitting schemes associated with the operator 

nh h 

to which we associate the numerical solution 

(n,,,oc^^^ov,^)"(zO). (3.14) 

Obviously, for r] = 0, II^^s is the identity operator. 

In the following, we show a normal form result on the abstract splitting 
method 

fHo ° "PP 

and then draw some dynamical consequences for the discrete solution (3.14). 



4 Statement of the result and applications 
4.1 Main result 

Let j = {ji, . . . ,jr) G , and denote by jj = {ai,Si) £ AAx {±1} for z = 1, . . . , r. 
We set 

= 6iUJai H h 5rUJar- 

We say that j = (ji, . . . ,jr) £ Ir depends only of the action and we write j G Ar 
if r is even and if we can write 

Vi = l,...r/2, ji = {ai,l), and ji+r/2 = {ai, -I) 

for some E M. Note that in this situation, 

— lai ■ ■ ' -^0^/2 

where for all a e M, 

Ia{z) = CaVa 

denote the action associated with the index a. Note that if z satisfies the condi- 
tion zj = 'Zj for all j G Z, then we have laiz) = |CaP- For odd r, Ar is the empty 
set. 

We will assume now that the step size h satisfies the following property: 
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Hypothesis 4.1 For all r G N, there exist constants 7* and a* such that \/ N £ 
N* andyj = (ji,...,jV) ^ A, 

|ji|,...,bV|<A^ =^ |l-e^''^(^-)|>|^. (4.1) 

Theorem 4.2 Assume that P £ T and h < Hq satisfies the condition (4.1). Let 
r G N* 6e /zxed. T/ien there exists a constant sq depending on r such that for 
all s > sq, there exist constants C and Eq depending on r and s such that the 
following holds: For all e < Eq and for all G V2s i^go,1 such that Hfj^gZ^ = 
with r] = e'^+^Z^ and 

< e and |b°|L < 1, 

II II ^ — II II 2s — ' 

if we define 

= (n^,, o (^^^ o <^^)"(zO) with r? = e'-+V4 (4.2) 
then we have is still real, and moreover 

||z"||^<2e for ri<-^, (4.3) 

and 

5^|a|2^|I,(z'^)-/,(z°)| <e^/2 for n < (4.4) 

aeJ\f 

The proof is postponed to section 5.3. 

Remark 4.3 As r is arbitrary, the condition rj = e'^+^Z'* implies that II^^s is 
gr+1/4 f^iggg j-Q ^/jg identity in Vs (cf (3.13)J. From the practical point of view, 
we can always assume that is beyond the round- off error, so that we can con- 
sider that (4.2) coincides with the numerical solution associated with the splitting 
method. The full understanding of the real numerical phenomenon taking into 
account the round- off error is clearly out of the scope of this paper. 

Remark 4.4 The condition \\zo\\^^ < 1 together with IT^^sz" = z^ implies that 
z^ = for j large enough which is actually the assumption we need. 



4.2 Verification of the non resonance condition 

We assume that the frequencies uJa, a ^ fulfill the following condition (in the 
next section we will verify this condition in different concrete case): 

Hypothesis 4.5 For all r G N, there exist constants "f{r) and a{r) such that 
ViVGN* andVj = (ji,...,j,) ^ A, 

bi|,...,|jV| <iV ^\^{j)\>J-. (4.5) 
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The next result shows that under the previous hypothesis, the condition (4.1) 
is generic. See [17, 16] for similar statements. 

Lemma 4.6 Assume that Hypothesis 4-5 holds , and let Hq and r be given num- 
bers. Let 7 and a be such that (4.5) holds and assume that 7* < (2/7r)7, 
a* > a + ma + r with a > 1 and m the constant appearing in (3.10), then 
we have 

7 -, . 

meas{ h < h^lh does not satisfy (4.1) } < C — 

7 

where C depends on a and r. As a consequence the set 

Z^Hq) = { /i < /iQ I /i satisfies Hypothesis 4-^} 
is a dense open subset o/(0,/io)- 

Proof. Denote 

R(/io, 7*, a*) = {h < holh does not satisfy (4.1)}. 
Assume that h G R(/io, 7*, «*)• There exist > 1 and j ^ Ar such that 

lii I , • • • , I jV I < and 1 1 - e*^^(^') | < 



For this j, there exist an £ G Z such that 

|1 _ gihnU)^ > l\27r£ - hn{j)\. (4.6) 
vr 

If ^ = 0, the previous inequality and (4.5) imply 

' ' ~ vr A^" 

which is impossible with the assumptions on 7* and a*. Hence, we can assume 
i ^ 0. Eqn. (4.6) implies 



m3)\ 



IT 



lire , 
h 



hy 

< 



and using (4.5) 

2Tri , , Httj* 1 



h\ < 



n{j) ' - 27 AT"*-"' 
Moreover, we have for this i 

\2Tri-hn{j)\ < TT 

whence using (3.10) 

2tt\£\ < vr + ChoN"^ 
where C is a constant depending on r. This implies 

W - i < i'.oA"". 
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Hence, R(/io,7*,a*) is included in the union of balls of center 

2TTf 

with |jl|,...,|ir| < iV, \l\<TT + ChQN'^, £y^0 
hoTX^* 1 

Hence, we have for a > 1 

/iovr7* 1 



and radius 



meas 



(R(/io,7*,«*))< E E 



27 AT"*-" 

\'j^\<N \t\-l<C.h^Nm 

^ ^ ^ / 1 \'^/iovr7* 1 (Ch^Y 



\h\<n 



iy 27 iV"*-Q-m'T V 27r 



— 7 jy^ 



Furthermore 



meas(n^*>oR(/iO) 7*) a*)) = 
and thus Z{h{)) has full measure. ■ 

4.3 Examples 

In this section we present two examples, other examples like the Klein Gordon 
equation on the sphere (in the spirit of [3] ) or the nonlinear Schrodinger operator 
with harmonic potential (in the spirit of [13]) could also be considered with these 
technics but will require some additional efforts. 

4.3.1 Schrodinger equation on the torus 

We first consider non linear Schrodinger equations of the form 

idtilj = -^ij + Vi^i) + d2g{ip,ij), rrET'^ (4.7) 

where V G C'^{T'^, M), g G C^iU, C) where U is a neighborhood of the origin in 
C^. We assume that g{z,z) G M, and that g{z,z) = 0{\z\^). The corresponding 
hamiltonian functional is given by 



d 



Let 4'a{x) = e*"'^', a G Z'^ be the Fourier basis on L^(T'^). With the notation 

/ 1 Nd/2 _ / 1 Nd/2 

^ = [^) 1^ CaMx) and Tp = l^—j 2^ rjacpaix) 
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the hamiltonian associated with the equation (4.7) can be (formally) written 



V) = ^aCa??a + ^ ^ i^ab Cai • • • Cap7?fei • • • ^b, • (4.8) 

aeZ'i r>3 a,b 

Here coa = + Va satisfying (3.10) with m = 2 and are the eigenvalues of the 
operator 

■0 1— > — Alp + y * 

Note that in (4.8) the sum is made over the set of multi-indices 

{{a,b) = {ai,...,ap,h,...,bq) e X {Z'^y with p + q = r 

and ai + • • • + Up — bi — . . . — bq = 0} , 

which corresponds to the set (3.6) in variables {zj)j^z (here we set Af = TH^). 

The relation -ff (^, ^) € M is equivalent to the fact that the coefficients P^h 
satisfy P^h = Pha which correspond to the hypothesis (3.8) in variables {zj)j^z. 
The fact that the nonlinearity P belongs to T can be verified using the regularity 
of g and the properties of the basis functions (pa, see [12, 4]. In this situation, 
it can be shown that the Hypothesis 4.1 is fulfilled for a large set of potential V 
(see [4] or [12]). 

The numerical implementation of splitting method is very easy in the case of 
Eqn. (4.7): The part corresponding to the equation 

idtip = —Alp + V -kip 

is easily solved in terms of Fourier coefficients, while the non linear part 

idttp = d2g{ip, tp) 

is a simple differential equation with fixed x G T'^. The use of a fast Fourier 
transform allows to compute alternatively the solution of the linear part and the 
solution of the non- linear part. 

4.3.2 Wave equation on the circle 

We consider the wave equation on the circle 

Uu-Uxx+mu = g{u), x G T\ t G M, 

where m is a non negative real constant and g a smooth real valued function. 
Introducing the variable v = ut, the corresponding hamiltonian can be written 

H{u, v) = I ]-{v^ + ul + mu^) + G{u) dx, 
15 



where G is such that duG = g. Let A := {—dxx + m-)^/^, and define the variables 
{P,q) by 

q ■= J^l'^u, and p = A'^l'^v. 
Then the Hamiltonian can be written 

Let lOa = -\/|ap + m, a G N =: be the eigenvalues of the operator A, and 
the associated eigenfunctions. Plugging the decompositions 

9(3^) = '^(la4>a{x) and = y^pa(/>Q(x) 

aGN aGN 

into the hamiltonian functional, we see that it takes the form 

aGN 

where P is a function of the variables -pa and qa- Using the complex coordinates 

ia = -^iSla + Wa) and r]a = ^{qa - iPa) 

the hamiltonian function can be written under the form (3.9) with P £ T (see 
[2, 12]). As in the previous case, it can be shown that the condition (4.5) is 
fulfilled for a 1 set of constant m of full measure (see [4, 2]). 
In this situation, the symmetric Strang splitting scheme 

h/2 h h/2 

corresponds to the Deuflhard's method [15]. Considering now the Hamiltonian 

H{z)=Ho{z)+PmhHo)z) 

where ^(x) a smooth function that is real, bounded, even and such that ^*(0) = 1, 
and where 



then the splitting schemes associated with this decomposition coincide with the 
(abstract) symplectic mollified impulse methods (see [16, Chap. XIII] and [9]). 
The fact that <I> is bounded makes that the functional z 1— > P(^{Hq)z) obviously 
belongs to T. 
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5 A normal form result 



5.1 Normal form 

Definition 5.1 Let N > be a real number. For a given multi-index j G Z'^ , 
let ip be the permutation such that 

\jh \ < ■■■ < 

We define the set 

Jr{N) = {j G Ir liij < (r - l)iV and < N}. 

Lemma 5.2 Let r > 3, and assume that j ^ JriN). Then j contains at least 
two indices with modulus greater than N. 

Proof. Let j = {ji, ...,>) G lr\Jr{N). We have M{j) = where M{j) is 
defined in (3.5). Assume that there exists only one index of modulus greater than 
A^. We can assume that | ji | > N and hence all the other indices are of modulus 
< (in particular, with the previous notation, we have ji = ji^). Hence we 
have 

|ji| < U2| + --- + |>| < (r-l)Ar 
and this implies that j £ J'r{N) which is a contradiction. ■ 

We motivate now the definition of normal form terms we introduce in the 
sequel. For a given number A^ and z G "P^ we define 

Nf (Z) = l«l'^ear?a 

\a\<N 

and 

|a|>Ar 

SO that 

||z||; = Nf(z) + Rf(z). 



Proposition 5.3 Let A^ G N and r > 3. Assume that the homogeneous polyno- 
mial 

jeIr\Jr{N) 

defines an element oJTJ^^''^ for some constants M and u. Then we have for all 
s > 2iy + i, M > s + 2 and for all z G VsiC), 

|{Nf ,Z}(z)| < Co\Z\^m,.N''~^^+''/'-^\\z\\'''-'r^{z). (5.1) 
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and 

Va G AT, \a\< N, \{Ia, Z}\ < Cq\Z\^m,. N^+'^-'\\z\\''~^ Rf (z). (5.2) 
Moreover 

|{Rf , Z]{z)\ < Co\Z\^M,. Ikli;-' Rf{z) (5.3) 
where Cq is a constant depending on s, r and the dimension d of M = N'^ ofTL'^. 
The proof of this proposition is given in the Appendix. 

Definition 5.4 An element Z G T^ ^'^ is said to be in normal form if we can 
write it 

r 

^=3 {3eAeUlc\Je{N)} 

In other words, a normal form term either depends only on the actions or 
contains at least two terms with index greater than N (cf. Lemma 5.2). 

5.2 Statement of the normal form result 

In the following, we set 

Bs{p) = {zers\\\z\\^<p}. 

Theorem 5.5 Assume that P £ T and h < ho satisfies the Hypothesis 4-1- 
Let ro > 3 be fixed. Then there exist constants sq, (3 and Nq such that for all 
s > sq, there exists constants C and for all N > Nq there exists a canonical 
transformation t from Bs{p) into Bs{2p) with p = {CN)~^ satisfying for all 
z £ Bs{p), 

\\t{z) - z\\^<{CNf\\z\\l and ||r-^(z) - z||^ < (CiV)'^||z||^ (5.4) 
and such that the restriction of t to the high modes is the identity, i.e. 

{t{z)), = Zj for \j\ > (ro - l)iV. (5.5) 
Moreover, r puts (p^ in normal form up to order tq in the sense that 

where is the solution at time A = 1 of a non- autonomous hamiltonian hZ{X) + 
R{X) with 
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• Z{X) G C([0, 1], T^Q^'"^^) for some Mi and vi depending on P, ro, s and ho, 
and for all A G [0, 1], Z(X) is a real polynomial of degree r under normal 
form such that 

\Z{X)Lm,.., < {CN)f'. (5.7) 

• R{X) G C([0, 1], with p < {CN)-I^ has a zero of order ro + l at 
the origin and satisfies and for all z G Bs[p), 

VAg[0,1], <{CNf\\zf^. (5.8) 

The proof is postponed to section 5.4 and 5.6. We first verify that this normal 
form theorem has the dynamical consequences announced in Theorem 4.2. 



5.3 Proof of the main Theorem 4.2 

We now give the proof of Theorem 4.2. 

First, let us note that as the Hamiltonian functions Hq, P, Z and R are real 
hamiltonians, and by definition of II^^s (which is symmetric in ^ and 77), it is 
clear that there exist 1^" G C-'^ such that for all n, we have = (C"i'^")- 

Let ro = r, and let C = C(r, s), sq = so{r) and /3 = /3(r) be the constants 
appearing in Theorem 5.5. We can always assume that 

so(r) > 2(r + l)/3(r). (5.9) 

Let s > so{r). Let Eq be such that 

< P.IO) 

Let e < eo- We define N such that 

(CAT)/? = £-1/2, (5,n) 

Notice that the assumption < e implies z^ G Bs{p) with p = {CN)~^. 

Furthermore we have II^^H^^ < 1 and together with H^^gz" = z^, this hypothesis 
implies that Zj=0 for j large enough. Actually let j G 2^ be such that |j| > 
(r — 1)A'", we have 

\f\\z]\<\jr<{{r-l)N)-'. 
i_ 

using (5.11) we have N = C~ e ^/s and hence 

< {r-l)-'C'e^. 

Now condition (5.9) implies that ^ > r + 1 and hence (as we can always assume 
that e < 1), 
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Therefore we get using (5.10) 

l/ll^l < e'"+^/'((r - < £-+1/2. 

As Uri^sZ^ = and 77 = e''+i/^, this imphes that 

V|j| > (r-l)Ar, 20 = 0. 

Let r defined by Theorem 5.5, and let = t-1(z"'). As r is the flow of a real 
hamltonian, there exist G such that y" = (C", (^") for all n. By definition, 
we have 

Vn>0, 2/'^+i = (r-ion^,,or)o((^^^oV)(y"). (5.12) 

and as r is the identify for high modes (see (5.5)), we have = for \j\ > 
(r - 1)A^. 

Using the definition of A^, the transformation r in the previous Theorem 
satisfies (taking p := 2e < -v/e): for all z such that \\z\\ < 2e, 



\t-^z)-z\\ < £-1/2 
I \ / " s ~ 



\Z\ 



|2 



< 4£3/2 (5.13) 



< b 

provided £0 is sufficiently small. Hence, we have = ||t-1(z'')|| < |£. 

We will show by induction that the following holds for all n E N: 

(i) < ||?/°||J + 2n£'^+i/8 

(ii) = for j>{r- 1)N. 

These assumptions are satisfied for n = 0. Assume that they hold for n > 0. 

Let ip be the application defined by Theorem 5.5 and ^''^(•2) be the flow 
associated with the Hamiltonian hZ{X) + R{X) defining the application -0 for 
A = 1. 

Using the results of Lemmas 5.6 and 5.7 below, we easily see that there exists 
a constant c depending on r such that for all A G [0, 1], ||V''^(y"')ll^ ^ c£. 

Let Ni = {r — 1)N. We have by hypothesis that R^^{y"') = 0. Furthermore 

Thus using the equation (5.3), (5.7) and (5.8) we get 

ARfl(V,^(y-)) 
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for some constant Ci depending on r and s. Hence we have 
d 



dA 



where Ci depends on r and s. Using the Gronwall Lemma, we obtain for all 
A G [0, 1] 

1/4 1/2^ 

We can always assume that Eq Cie^o < 1. Hence we get 

VAg[0,1], RfHV'^(y")) < e"^^/^- 

On the other hand, using (5.1) we have 
d 



< 



CiiV^+'^+2+'^/2-^Rfi(V;^(y"))(||V;^(z) 



s 



for some constant Ci depending on r and s. We can always assume that s > 
j3 + V + d/2 + 2. Using the previous estimates, we get 



dA 



NfnV'"(y")) 



and then 



Nfnv^'(y"))<l|y"ll! + Cie'^+'/'- 



Let = 9?^^^ o As for all 2; we have ||z||^ = N^i(z) + R^i(z) and as (^^^^ 



preserves all the actions, therefore 



\y 



(5.14) 



Now by construction (cf. (5.12)) 

y"+'=r-ion,,,or(y"). 
As T is the identity for modes |j| > A^i, we have 

Rfn^(y")) = RfHy") < e"^^^^ = ^• 

Hence by definition of the projection H^^g we get that 

(n^,, o Tir))^ = 0, for \j\ > (r - l)N = N,. 
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As r ^ is the identity for modes greater than (r — 1)A^, this shows (ii) for n + 1, 
i.e. 

y^+^ = for \j\ > (r - l)iV. 
Let z be such that Zj = for \j\ > (r — 1)A^. We have 

lil<Af 

since we can always assume /? > id. 
Writing 

o n^^^ o r = I + o (n.;y^s — /) o T 
and as r leaves the set (zj)|j|>jv^ invariant, we get using (5.4), 



Thus we get 



||yn+l||2 ^ ||^„||2 ^ ^^,+1/4 ^ 3^.+l/8 < ||y„||2 ^ 2^,+i/8_ 

This shows (i) for n + 1. 

In particular, for all n < e"^"*"^ we have (recall Hy^H^, < |e) 

||y"||;<(|e)' + 2s2+V8 
and hence (provided Eq is small enough) 

Now using (5.13) for the application r, we easily see that \\z"'\\ < 2e as long as 
n < e*""^. This proves (4.3). 

The proof of (4.4) is obtained similarly using (5.2) and we do not give the 
details here. 



5.4 Formal equations 

We give now the general strategy of the proof of the normal form Theorem 5.5, 
showing in particular the need of working with non autonomous Hamiltonians 
and of considering the non resonance condition (4.1). 

We consider a fixed step size h satisfying (4.1). We consider the propagator 

V>Ho °Vp = ^Ho ° fhP- 
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We embed this application into the family of applications 

^Ho°'PhP, AG [0,1]. 

Formally, we would like to find a real hamiltonian x = x(A) and a real hamilto- 
nian under normal form Z = Z{X) and such that 

VAG[0,1] v'h,o^^^pOip^^^^=^^^^^oip%^oip>^z^^y (5.15) 

Let G T'si'C) and z'^ = Lp^^^^ o c^^p o 99^^^^ (z*^). Deriving the previous equation 
with respect to A yields 



Using Lemma 3.2 that remains obviously valid for non autonomous hamiltonian, 
we thus have 



dz^ 

where ^(A) it the time dependent real hamiltonian given by 

A{x) = hP o (^^j + x{x) o (^^^ o 

Using the same calculations for the right-hand side, (5.15) is formally equivalent 
to the following equation (up to an integration constant) 

V A G [0, 1] /iP o (^-J + X o 99-^ o (^-J = x(A) + hZ{X) o ^-^^^ o (^-J. (5.16) 

which is equivalent to 

VAg[0,1] xWoip'k,-xWoip^^ = hP-hZ{X)o^-^^y (5.17) 

In the following, we will solve this equation in x(A) and Z{X) with a remainder 
term of order r in z. So instead of (5.17), we will solve in section 5.6 the equation 

VAg[0,1] x(A)ov9^„-x(A)o(^-^ = /iP-(/jZ(A) + i?(A))o(^--\^. (5.18) 

where the unknown are x(A), and Z{X) are polynomials of order r, with Z under 
normal form, and where R{X) possesses a zero of order r + 1 at the origin. 
In the following, we formally write 

r r 

X(A) = J2 ^m(^) -=5^5^ 
e=3 £=3 jell 
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and 

r r 
^=3 £=3 j&li 

where here the coefficients Pj{X) are unknown and where X[e\{^) s-^icl Z[^](A) 
denote the homogeneous polynomials of degree ^ in x{^) ^-iid Z[\). 
Using the assumptions on P, we can write 

r 

P = A + B = Y,P[t]+B 

where A E T^''^ and B £ 7i^{Bs{po)) for s > sq and po sufficiently small. 
Moreover, B has a zero of order r + 1 at the origin. 

Identifying the coefficients of degree £ < r in equation (5.18), we obtain 

Xm(A) o LfH^, - Xm(A) = hP[i] - /i%](A) + /iG[f](A;x*,-P*,^*)- 

where G is a real hamiltonian homogeneous of degree I depending on the poly- 
nomials X[fc]) P[k] s-nd Z\^^ for k < I. In particular, its coefficients are polynomial 
of order < ^ of the coefficients Xji Pj a-^d Zj for j £ I/^, k < £. 
Writing down the coefficients, this equation is equivalent to 

yjeir (e^'^^W) - l)xj = hPj - hZj + hGj 

and hence we see that the key is to control the small divisors e^^^^^'^ — 1. 



5.5 Non autonomous Hamiltonians 

Before giving the proof of Theorem 4.2, we give easy results on the flow of non 
autonomous Hamiltonian. Let Q{X) G C([0, 1],7^^^''^) for some r > 3, M > and 
1/ > 0. We set 

IIQIIt-a/.^ = max \Q\^M.u . 
ag[o,i] 

The following results extend the properties already proved in [12] or [2] and 
needed in the proofs below. 

Lemma 5.6 Let /c G N, M G N, G [0,oo), s G M with s > + 3/2, and let 
P(A) G C([0, l],Ti^^'^) be a homogeneous polynomial of order k + 1 depending on 
A G [0,1]. Then 

(i) P extends as a continuous polynomial onVs{C) depending continuously on 
A G [0, 1], and there exists a constant C such that for all z G 'Ps(C) and all 
AG [0, 1], 

|P(A,z)| <C||P||^M,.||z||'+'. 
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(ii) Assume moreover that M > s + 1, then the Hamiltonian vector field Xp(^x) 
extends as a bounded function from Vs{C) to Vs{C) depending continuously 
on \ ^ [0,1]- Furthermore, for any s > v + \, there exists a constant C 
such that for any z G Ps(C) and A G [0, 1], 



||Xp(A)(^)ll,<C||P||^M..||z||'. 

Lemma 5.7 Let r > 3, M > and let 

r 

Q(A,z) = ^^Q,-(A)z,- 

£=3 jeii 

be an element of €{[0,1],%^''^). Let ^q(^x) A^'^ associated with the non 

autonomous real hamiltonian Q{X). Then for s > i' + 3/2 there exist a constant 
Cr depending on r such that 

p<mf{l/2,Cr\\Q\\~i,.) =^ VAg[0,1], ^^^x)iBs{p)) C Bs{2p). 

(5.19) 

Moreover, if F{\) G C([0, l\^'H'^{Bsi2p)) has a zero of order r at the origin, then 
F{X) o (Pq(^x^ has a zero of order r at the origin in Bs{p). 

Proof. Let z'^ = (^^^^^(a:'^). Using the estimates of the previous lemma, we 
have 

< Cr\\Q\\^M,A\A\,{\\zX + \\'%~') 

for some constant Cr depending on r. Hence, as long as < 1, we have 

^\\zX<'^Cr\\Q\\^MA\A\l- 

By using a comparison argument, we easily get that for z^ G Bs{p) we have 

VAg[0,1], ||z^||^<2||z0||^. 
This shows (5.19) and the rest follows. ■ 

The next result is a consequence of Prop 6.3 in [12]. The only specificity is 
the control of the sum of the indices, and the evolution of the norm 

Proposition 5.8 Let ki and k2 two fixed integers. Let P and Q two homoge- 
neous polynomials of degree ki + I and k2 + I such that P G C([0, 1],T^^^) and 

Q G C([0, l],'7;f+7) for some ui > 0, 1^2 > and M > 0. 
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Then {P, Q} defines a homogeneous polynomial of degree ki +k2, and for all 
M' and v' such that 

M' < M - max(z^i, ^2) - 1 and v' > U1 + U2 + 1, 

we have {P, Q} G C([0, 1], 7^^^'^_^ ) and 

\\{P^Q}\\tM'.m' < C\\P\\^M.ui ||(5||^A/,^2 

for some constant C depending on M , v , M' , v' , ki and /c2- 

Proof. The proof is clear using Proposition 6.3 in [12]. We only need to 
verify the fact that the summations are always made over set of indices with zero 
moment which is trivial. ■ 

Lemma 5.9 Let xW an element of C {[0,1], tJ^^''^) for some M > and 
V > 0. Let t(A) := '^^{^x) /iloiy associated with the non autonomous real 

hamiltonian x(A). Let g G C([0, 1],%^^''^), then we can write for all ctq G [0, 1], 

g{ao) o r(cro) = g{cro) 
+ /Z ••• / Ad^(<,^.)0...oAd^(^^)^(c7o) dcTi---dfTfc + i?(ao) (5.20) 
where by definition Adp{Q) = {Q,P} 
Ricro) = J ■ J (Ad^(^^) o • • • o Ad^(^j)C/(cro)) o r(o-r)d(Ti • • -do-,.. (5.21) 

Each term in the sum in (5.20) belongs (at least) to the space C([0, 1],?^^,^ ''^ ) 
where 

J/' = (r + l)(i/ + 2) and M' = M - v' . 

The term R{cro) defines an element of 7i'^[Bs{p)) for s > v' + 3/2 and p < 
inf(l/2, (7^11x11^-1/,^) and has a zero of order at least r + 1 at the origin. 

Proof. For a fixed (Tq G [0, 1], we have 

^5(0-0) o t(A) = {(/(fJo), x(A)} o t(A). 

Hence, we have that 

3(0-0) o r(cro) = g{ao) + {kd^^^„^)g{aQ)^ o r(cri) dai. 
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Repeating again the same argument, we have 

g{ao) o t(cjo) = giao) + / (Ad^(^^)5f(o-o))dcJi+ 

Jo 



Jo L (^^^('^2) ° Ad^K)5'(^o)) oT(o-2)do-idc72 



The equation (5.20) is then easily shown by induction. The result then follows 
from the previous propositions. ■ 

For a given polynomial x G C([0, 1],%^^'^) with r > 3, we use the following 
notation 

r r 

X{\z) = Y.X[^]W = Y.Y. (5.22) 

where X[t]W 6 C([0, 1],%^^''^) is a homogeneous polynomial of degree I. 

Proposition 5.10 Letx{X) he an element o/C([0, 1],7^*^''^) for some M > anc? 
1/ > 0. Xe^ '/'^(•;)^) ^/je flow associated with the non autonomous real hamiltonian 

x(A). Let g G C([0, 1],tJ^^''^), then we can write for all X G [0, 1], 

5(A)o^^(,)=5M(A)+tM(A) 

where 

• ^^(A) e C~([0, 1],%'^^'"') with ui = {r + l)(i/ + 2) and Afi = M - z^i. 

Moreover, if we write 

r 

e=3 

where 5'[^](A) is a homogeneous polynomial of degree I, then we have for all 
£ = 3,...,r, 

S[i]{X) = g[i]{\) + Gye^{\;x*,g*) 

where G^p^{x*,g*) is a homogeneous polynomial depending on A and the 
coefficients Sj are polynomials of order < I of the coefficients appearing in 
the decomposition of g and x- Moreover, we have 



m=3 m=3 

(5.23) 

where C depends on i, M and v. 
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• rW(A) G n'{Bs{p)) for s > I/' + 3/2 and p < inf(l/2, a||x|LM,.) and 
has a zero of order at least r + I in the origin. Moreover, we have for all 
z G Bs{p), 

vag [0,1], i|x^w(A)(^)IL<a(x*,5*)lkn; 

where 

r r 

\\9im]\\j-M,-^ ) (i + ^ IIxhIIt*^'" ) 

m=3 m=3 

with C depending on r, M and v. 

Proof. Using the previous lemma, we define S'^^'^ as the polynomial part of 
degree less than r in the expression (5.20). The remainder terms, together with 
the term -R(A) in (5.21), define the term T('"^(A). The properties of S^'''^{X) and 
T^''^(A) are then easily shown. ■ 



5.6 Proof of the normal form result 

Proposition 5.11 Let P G T and N be a fixed integer. Let M > u' := (r + 

+ 2) and M' = M - v' . Then there exist 

• a polynomial x £ C([0, 1],7^^^ ) 

x(A) = Exm(A):=^ XjWzj 

• a polynomial x £ C([0, 1],%^^ ) 

r r 

Z(A) = j;Z[,](A):=5^ ^i(A)^i 

^=3 1=3 {j(iAiVJli\Ji{N)} 

under normal form, item a function R{X) G C{[0,l],Ti.s{Bs{p))) with p < 
cqN~^ for some constant cq > and [3 > 1 depending on r, M, P, and 
having a zero of order r + 1 at the origin 

such that the following equation holds: 

VAg[0,1] xWov'k,-xWo^^^ = hP-ihZ{X)+RiX))o^-^^^. (5.24) 
Furthermore there exists Cq > depending on P, v , r, M such that 

\X\q-M' y + \Z\ 

,j-M' ,v' 
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and moreover 

VAg [0,1], \\Xnix){z)\\^<CoN^\\z\\l 
for z G Bs{p) with p < cqN~I^ . 

Proof. Identifying the coefficients of degree £ < r in the equation (5.24), we 
get 

X[i] ° Vho - Xl£] = hPii] - hZ[f,^ + hGif^ix*, P*,Z^). 

where G is a real hamiltonian homogeneous of degree £ depending on the poly- 
nomials X[k]i P[k] s-rid ^[k] foi" k < i. In particular, its coefficients are polynomial 
of order < £ of the coefficients Xj, Pj and Zj for jf S Xfc, fc < £ and satisfy- 
ing estimates of the form (5.23). Writing down the coefficients, this equation is 
equivalent to 

yj€lr {e'^^'^^) _ i)^. = hp. _ hZj + hGj. 
We solve this equation by setting 

Zj = Pj + Gj and Xj = for j ^ Ai>VJ Ii\Ji,{N) 

and 

Zj = and = -^^jA_(p. + G,.) for j e Je{N)\A,. 

Using (4.1) and the result of Proposition 5.10 we get the claimed bound for some 
/3 depending on r. 

To define R, we simply define it by the equation (5.18). By construction, it will 
satisfies the announced properties. ■ 

Proof of Theorem 5.5. Integrating the equation (5.18) in A, it is clear that 
the following equation holds: 

V A G [0, 1] o o (^J(^) = ^J(^) o ip%^ o ip^z(\)+R{X) ■ 

Note that using Proposition 5.10 and (5.19) we show that for s > z/' -|- 1 and 
z G Bs{p) with p = cN~^ we have 

K^>.){z)-zl^<CN^\\z\\l. 
This implies in particular that 

ll-IL^II^J(A)(-)IL + ^^-'11-11. 

For sufficiently large, this shows that '^'^(^x^ is invertible and send Bs{p) to 
Bs{2p). Moreover, we have the estimate, for all A G [0, 1], 

\\Kw)~\z)-z\\^<CN^\\z\\]. 

We then define r = V'^(a) = ^\z{\)+R{\) ^^'^ verify that these application 

satisfy the condition of the theorem. ■ 
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6 Appendix: Proof of Proposition 5.3 



Let j G IrX'JriN). It is clear that for a G M, we have, with the notation 
unless (a, 1) or (a, —1) appears in j. Moreover, if this is the case, we have 

\{Ia,Zj}\<2\zj\. 

where we set 

l^jl ~ Nil I ■ ■ ■ I^Jrl 

for j = {ji, ■ ■ ■ , jr) G ■ Hence we can write 

|{Nf,Z}(z)|<2 \Zj\\zj\ 

\k\<N {jeIr\JriN)\jDk} 

where k = (a, ±1) G Z in the first sum. We thus get using (3.7) 



|{Nf,Z}(z)|<2|Z|^„,. Y 



1 2s/ 



Zi ■ 



\k\<N {j£lr\JriN)\jDk} 

Using Lemma 5.2, the indices in the previous sum are such that at least two of 
them are greater than N. As |A;| < A'^, these indices cannot be equal to k. Hence 
we can rewrite each indice j containing k as {j,k) where j G Z^~^ contains 
at least two indices greater than N. Using the symmetries in the sum, we can 
moreover assume that the indices are ordered in such a way that \ji\ > |j2| > . . .. 
Hence, we can rewrite the previous sum as 

|{Nf,Z}(z)|<C|Z|^M.. Y ^^-^^ 

jeZr^l,\jl\,\j2\>N,\k\<N ' 

where fi{j, k) and S{j, k) denote the values of /i and S associated with the r-uple 
(ji) j2> ■ ■ ■ ,jr-i, k). Here, C denote a constant depending on r. 
Since k) < S{j, k) and /u(j, k) < |j2| we have for M > 2 



\{N^,Z}iz)\<C\Z\ 



% 



X 

iG2r-i,bi|,|j2|>Af, |fc|<7V 

Then use \k\ < \ji\ to obtain 



\{K,z]{z)\<c\z\ 



X ^ 



j(iZr-l,\h\,\j2\>N,\k\<N 
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By Cauchy Schwarz, one has for s > 1/2 

Ei^'i<ii^ii^fEi^r'') (6-3) 

and thus we get from (6.2) 

\{N^,Z}{z)\<C\Z\^M..\\z\r^ 

^ 2 

X 



\ji\,\j2\>N, \k\<N 

Hence, introducing the sequence {bj = \ j\''\zj\)j^z £ ^^(-^) we can write 
\{N^,Z}{z)\<C\Z\^.,.\\z\r' 

^ 2 

X 



\ji\,\j2\>N,\k\<N 

Moreover, the sum in < in (6.4) yields by Cauchy-Schwartz inequahty 



^ fefc < CiV^/y N^(z) < CN'^/^W 

\k\<N 



Z\\ 



where d is the dimension of = N'^ or Z"^. Hence, we get from (6.4) using 

|j2|>iV 

|{Nf ,Z}(z)| < CiV-+^+^+'^/2|Z|^M,. \\z\r' J2 ( i , |-|_|- 

Iii|>li2|>iv ^ ^l-^il 1-^21/ 

and this concludes the proof of (5.1) since, if a and c are two sequences in i'^{2) 
we have by a convolution argument 

+ - \l\ ) l''^'' - MiHz) (6-5) 

for some universal constant C. 

Note that (5.2) is an easily shown by similar calculations (the only difference 
lies in the fact that there is no summation in < A^). 

We now show (5.3). 
As contains only indices greater than N, we can write (see (6.1)) 

|{Rf,Z}(z)|<C|Z|^„,. J2 IfcP^ g/- Lm l^jll^'^l (6-6) 

jeIr-i,\ji\>N,\k\>N ^ 
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where the sum is made over ordered indices \ji\ > |j2| > •••• Note that in 
opposition with the previous situation, we cannot ensure that \j2\ > in this 
sum. We first notice that, for all k and j, 

l*l|||<2b,|, (6.7) 

Actually, if \k\ < 2ji then (6.7) holds true since < 1- Now if > 2ji then 

S{l,j) > \\k\ - \ji\\ > and thus 

|A:|||^<2Mj,A:)<2Ui|. 

Then we distinguish two cases in this sum (6.6): 

|{Rf,Z}(^)| <C|Z|^M,. (/1 + /2) 

corresponding to the two cases \j2\ < |A;|, (/i) and \k\ < \j2\, {12)- 
Case 1: |j2| < |A;| 

In this situation, we use (6.7), fi{j,k) = \j2\, fJ-ijjk) < S{j,k) to conclude 
for M > s + 2 

oeIr-u\h\>N,\k\>N Vi + llJll MJ 

Then use (6.3) and the notation {bj = \ j\^\zj\)j^z G ^^(^) to get 

< C\\z\\lr^R^ {z) 

where we have used again (6.5) for ^^-^ ^ and (6.3) for Ylj2 ' 
Case 2: |j2| > \k\ 

In this situation, we still have iJ,{j,k) < |j2| and using that both \ji\ and \j2\ 
are greater than we get for M > 2 

h < C\\zf~' \ji\W^'\k\'^' , ,. |_|. X Ihl'^-'lzMzj^l 

\ji\>NAj2\>\k\>N ^ "^l-^^l 1-^21^ 



lii|>A^, \j2\>\k\>N 

<C\\z\CR^{zf/' 



i + lii|-U2l7 \k\^/^ \j2\'/^-^ 
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where in the last inequahty, we used that bj-^ , and y^^s/2-2-L^ respectively 

in f{Z), f{Z) (for s > 1) and f{Z) (for s > 4 + 2z/) and we used again (6.5) 
and (6.3). 
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